Polarization Enhancement in Short Period Superlattices via Interfacial Intermixing 
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The effect of intermixing at the interface of short period PbTiOa/SrTiOa superlattices is studied using first- 
principles density functional theory. The results indicate that interfacial intermixing significantly enhances the 
polarization within the superlattice. This enhancement is directly related to the off-centering of Pb and Sr cations 
and can be explained through a discussion of interacting dipoles. This picture should be general for a wide range 
of multicomponent superlattices and may have important consequences for the design of ferroelectric devices. 



Ferroelectric superlattices are a promising paradigm for 
creating novel materials for device applications. In ideal su- 
perlattices with perfectly flat, compositionally abrupt inter- 
faces, first-principles calculations and experimental studies 
have shown how factors such as strain due to lattice mis- 
matchea'.^.3iV.''.A^-9-io.i 1,12,13,14^ ^^^^^ compensation'^i^^ 

and bonding at the interfac&i^ can be controlled to enhance 
the ferroelectric properties of the superlattice. Work by Daw- 
ber and coworkers show that the behavior of very short period 
superlattices may deviate from that extrapolated from longer 
periods". While, these effects may be the result of a number 
of different factors, it seems likely that changes at the interface 
dominate at these length scales. Recent high resolution CO- 
BRA studies on SrTiOg (STO) supported PbTiOg (PTO) thin 
films further suggest that cation intermixing may be present at 
the interface of such superlattices^". 

Knowledge of how interfacial cation intermixing influences 
the aggregate properties of a superlattice can be useful in 
tuning its properties. As the superlattice period decreases, 
these effects become more important and, depending on their 
strength, may dominate. In this paper, we present a first prin- 
ciples study on the effect of interfacial cation intermixing on 
the polarization of short period PTO/STO superlattices. Our 
results demonstrate that intermixing significantly enhances 
the polarization and that this enhancement is directly linked 
to larger Pb displacements in the intermixed layers. 
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FIG. 1: Representative 2 PTO / 2 STO structures with intermixed 
interfaces. Grey squares are PTO unit cells and white squares are 
STO. Intermixed layers have a total composition of 50% PTO / 50% 
STO. The two intermixed layers account for 1 PTO and 1 STO layer. 



We use density functional theory (DFT) to examine the ef- 
fects of interlayer cation mixing on the polarization within 
PTO/STO superlattices. We compare superlattices with sharp 
interfaces with three different compositionally equivalent in- 
termixed systems of m=l-3 PTO layers and n=l-3 STO layers 
(figure [1]). Checkered, striped and 25%/75% intermixed lay- 
ers were studied. Intermixed layers were placed at each of the 
mixed PTO/STO interfaces; giving a total composition of 50% 
PTO / 50% STO. The two intermixed layers account for one 
PTO and one STO layer. All calculations were performed us- 
ing projector augmented wave (PAW) potentia b^^'^^ with the 
Vienna Ab initio Simulation Package (VASP v4.6.26)^', with 
the local density approximation for the exchange correlation 
functional. A 600 eV (22 Ha) cutoff and a 4x4x4 k point 
mesh were used. STO was found to be cubic with a lattice 
constant of 3.863 A (experiment: 3.905 A^"*). The tetrago- 
nal PTO lattice constants were computed as a=3. 867 A and 
c=4.033 A (experiment: a=3.904 A and c=4.152 A^^). This 
agreement is typical of LDA calculations for ferroelectric per- 
ovskites. In all calculations the in-plane lattice constant a of 
the superlattice was constrained to that of the STO substrate 
(3.863 A) while the c lattice vectors were optimized within the 
PAmm space group. AH ionic coordinates were fully relaxed 
until the Hellman-Feynman forces on the ions were less than 5 
meV/A. For bulk PTO constrained to the STO in-plane lattice 
constant of 3.863 A the c axis lattice parameter was 4.039 A 
and the polarization was 0.79 C/m^. Polarizations were com- 
puted using the Berry phase method^^. 

Figure |2] shows a comparison between the Berry phase po- 
larization and the polarization calculated using bulk Born ef- 
fective charges (Z*) and cation off-centering for all the in- 
termixed and unmixed superlattices studied. Cation displace- 
ments were computed from the center of their respective oxy- 
gen cages (i.e. 6 nearest oxygen neighbors for Ti and 12 for 
Pb and Sr). The bulk Z*s were: Zf,^ g^. = 2.7 and Z;^,; = 4.6. 
The remarkable agreement between both methods suggests 
that the bulk Z*s remain relatively unchanged in all super- 
lattices. This analysis allows us to relate changes in the total 
polarization to the off-center displacements of the individual 
cations. 

In figure[3]we see that the striped interface shows a marked 
increase in the total polarization as a function of percent PTO 
for nearly all short-period superlattices studied. While the 
striped interface is certainly the most dramatic example of 
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FIG. 2: A comparison between the polarization calculated using the 
Berry phase method and bulk Z*s and cation off-centering. The solid 
squares represent individual data points taken from the subset of cal- 
culations for PTO/STO superlattices with and without intermixing. 



polarization enhancement in an intermixed system, both the 
25/75 and checkered interfaces display significant increases 
in their polarization relative to the sharp interface, suggest- 
ing that this trend is the consequence of interfacial cation in- 
termixing. This polarization enhancement can be understood 
by considering the magnitude of the average cation displace- 
ments in both intermixed and sharp superlattices. Figure |4] 
depicts the average cation off-center displacements for Pb, Sr 
and Ti ions as a function of PTO concentration for superlat- 
tices with striped (right) and unmixed interfaces (left). These 
results indicate little change in the average Sr and Ti off- 
center displacements and clearly demonstrate that the changes 
in the polarization within the intermixed superlattices are due 
to changes in the off-centering of Pb cations. Once again, this 
trend is common to all intermixed configurations studied. 

To further elucidate the changes induced by interfacial 
cation intermixing we compared the changes in the average 
layer-by-layer cation off-centering for Pb and Sr cations in 
striped interfaces to those in pure layers (Figure |5]). Here 
we see that the magnitude of the average layer-by-layer Pb 
off-center displacements in unmixed Pb layers is essentially 
equal to that of Sr cations in unmixed Sr layers. In addition, 
the average layer-by-layer Sr cation off-center displacements 
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FIG. 3: Polarization as a function of PTO concentration for 
mpTo/nsTO short period superlattices. Solid squares (■) represent 
striped interfaces and open circles (o) are for sharp interfaces. The 
solid ( — ) and dashed lines(- -) are fits through the data. 



FIG. 4: Average cation displacements as a function of PTO concen- 
tration for the unmixed (left) and striped (right) interfaces. Squares 
(■) represent Pb data, circles (•) are for Sr and triangles (A) are Ti. 
Solid and dashed lines are fits to the data. Pb off-center displace- 
ments are greatly enhanced in superlattices with intermixing at the 
interface. 

decrease when comparing an unmixed layer to an intermixed 
layer Conversely, Pb cations off-center more in the inter- 
mixed layers. 

Now we turn to the discussion of how the differences be- 
tween unmixed and intermixed superlattices arise. First, we 
recall the present understanding of polarization in unmixed 
superlattices. PTO is a well-studied ferroelectric material with 
a double well potential which results in a spontaneous polar- 
ization in zero electric field. STO, on the other hand, is an in- 
cipient ferroelectric with a very shallow potential energy well 
centered around the paraelectric state. Previous first princi- 
ples superlattice studies have shown that in ideal ferroelectric- 
paraelectric superlattices with compositionally sharp inter- 
faces the ferroelectric layers, e.g. PTO, with polarization 
along the stacking direction, induce a polarization within the 
paraelectric layers. Furthermore, the magnitude of the layer- 
by-layer polarization throughout the superlattice was nearly 
constant2iJii2^; a direct consequence of the minimization of 
the electrostatic energy associated with the buildup of polar- 
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FIG. 5: Average layer-by-layer cation displacements as a function 
of PTO concentration in intermixed layers and unmixed layers for 
PTO/STO superlattices with striped interfaces. Solid squares (■) 
and circles (•) represent intermixed Pb and Sr layers. While, open 
squares (□) and circles (o) are for unmixed Pb and Sr layers. Solid 
( — ) and dashed lines (- -) are guides for intermixed and unmixed 
data, respectively. 
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TABLE I: AiJmixing relative to bulk for interfacial intermixed and 
unmixed superlattices. In all cases intermixing is preferred. 



FIG. 6: The percent enhancement in Pb cation off-centering as a 
function of number of Pb nearest neighbors. Shaded bars represent 
the average and standard deviations for superlattices with various pe- 
riods and compositions. The values in parentheses represent the per- 
cent PTO/STO in the layer The inset depicts the arrangement of the 
two types of Pb cations in the 75 PTO/25 PTO layer. Empty squares 
(□) are for Sr cations, filled circles (•) represent Pb cations with four 
nearest neighbors, and hashed circles are Pb cations with 6 in-plane 
neighbors (2 nearest neighbors and 4 second nearest neighbors). This 
plot was generated using data from all three intermixed superlattices 
as well as unmixed superlattices. 



ization charge at the interfaces^. The polarization of the super- 
lattice is then intermediate between zero and the full sponta- 
neous polarization of the ferroelectric layer material, increas- 
ing with increasing layer fraction of the ferroelectric phase. 
Our results for unmixed short-period PTO/STO superlattices 
are fully consistent with these previous results. In particular, 
the STO and PTO layers have nearly the same polarization, 
less than the full bulk spontaneous polarization of PTO, which 
increases with the compositional fraction of PTO. Since the 
Z*'& of Sr and Pb in STO and PTO are very similar, the Sr and 
Pb displacements are correspondingly nearly equal. 

Next, to understand the effects of intermixing we assign to 
each PTO or STO unit cell a single dipole, corresponding to 
the local polar distortion. Each layer is to be regarded as a 
square planar aiTay of dipoles. In an unmixed layer, where the 
dipoles are all PTO or all STO, the magnitudes of all dipoles 
in that layer are the same. On the other hand, in an intermixed 
layer where the potential energy wells governing the dipole 
moments of the two types of unit cells are significantly differ- 
ent, we expect the magnitudes of the dipoles to be different. 
Without changing the net polarization of the layer, the PTO 
unit cell dipoles can be larger, and the STO unit cell dipoles 
smaller, reducing the energy cost for maintaining the polariza- 
tion. As we discuss further below, depending on the details of 
the potential wells and the Pb/Sr arrangement in the layer, this 
gain in energy may drive an increase in the overall polariza- 
tion of the layer, and consequently of the superlattice. 

A more precise understanding of this argument can be ob- 
tained by examining the Pb off-center displacements as a 



function of the number of nearest in-plane Pb neighbors. To 
isolate the intermixing effect, we report these values relative to 
the average Pb displacements in compositionally identical un- 
mixed superlattices with the same period. As seen in Figure|6] 
the observed enhancement depends on the composition of the 
intermixed layer and the number of Pb in-plane neighbors in 
the layer First, there is a rough negative correlation between 
the Pb fraction in the intermixed layer and the enhancement of 
the Pb off-centering. Second, at a given Pb fraction there is a 
negative correlation between the number of Pb in-plane neigh- 
bors and enhancement of the Pb off-center displacements. 

These trends can be understood by considering the dipole- 
dipole repulsions within the intermixed layer Each pair of 
aligned Pb dipoles costs electrostatic energy, which increases 
as the spacing between the dipoles decreases or as the mag- 
nitude of the dipoles increases. A smaller Pb fraction in the 
layer decreases the relative importance of this contribution, 
favoring an increase in the magnitude of the Pb dipole. At a 
fixed fraction, the size of this energy will depend on the ar- 
rangement of the Pb dipoles in the layer, with a smaller num- 
ber of close (first and second) Pb-Pb neighbors decreasing the 
relative importance of this contribution and similarly favoring 
an increase in the magnitude of the Pb dipole. This final point 
is exemplified by considering the 75/25 mixed layer (see Fig- 
ure|6]inset). Here we see that there are two distinct types of Pb 
cations, the 4-Pb and the 6-Pb cation. Figure |6] shows that in 
this case the 4-Pb cation has a much larger enhancement, con- 
sistent with the decreases in dipole-dipole repulsions with de- 
creasing first and second Pb nearest in-plane neighbors. Sim- 
ilarly, Pb dipoles in striped (2-Pb nearest neighbors) layers 
are enhanced relative to those in the compositionally equiv- 
alent checkered (4-Pb neighbors) layers. In all cases, larger 
Pb dipoles in the intermixed layers increase the average po- 
larization of that layer, which, through the same electrostatic 
arguments presented for the sharp interface superlattices, fa- 
vors an enhancement of the total polarization. From this it can 
be understood that the striped interface superlattices, in which 
the Pb cations have fewer Pb neighbors than the checkered and 
the combined 25/75-75/25 superlattices, should exhibit more 
enhanced polarization, as observed in the first principles cal- 
culations. 
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Finally we consider whether interfacial cation intermixing 
is energetically favorable. TableUshows the Ai?mixing relative 
to bulk PTO and STO for the three interfacial intermixed cells 
and the sharp interface. In all cases, we see that intermixing 
is preferred over sharp interfaces. 

In conclusion, our results show that intermixed PTO/STO 
interfaces enhance the polarization of a superlattice relative to 
sharp interfaces. This enhanced polarization is strongly linked 
to larger Pb displacements in the intermixed layers and can 
be rationalized through a consideration of the dipole-dipole 



interactions. Since this model is based solely on the interac- 
tion of dipoles constrained by potential energy wells, these 
concepts should be general to all multicomponent superlat- 
tices in which there is a significant difference in the ferro- 
electric and/or dielectric properties of the substituent layers. 
This study, combined with knowledge from previous studies 
on strain and internal electric field effects in superlattices, sug- 
gests a new mechanism for tailoring short period superlattices 
for specific applications. 
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